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Abstract 

This  report  deals  with  the  heaving  and  pitching  motions  of 
unpropelled  ship  models  in  irregular  bow  or  stern  seas.  fine 
record  of  the  motions  is  obtained  from  a knowledge  of  the  /.uter 
surface  history  at  one  station  along  the  model,  using  a F-  rier 
integral  analysis,  as  oiscvsn.eu  by  Fuchs^2)*.  Thus,  time  nistories 
of  heaving  and  pitching  of  the  ship  model  can  be  expresses,  as  the 
convolution-type  integral  of  the  recorded  wave  motion  and  a kernel 
function.  This  kernel  is  the  Fourier  transform  of  the  ship  model's 
response  to  a sinusoidal  forcing  function.  Using  the.  Froude-Krileff 
hypothesis,  kernels  are  explicitly  computed  for  oscillations  of 
a rectangular  block.  The  kernel  for  a i;»uel  of  ship  form  is  de- 
termined from  the  results  of  a series  of  experiments  using  sinu- 
soidal waves,  and  independently  from  a single  experiment  using 
irregular  waves.  Predicted  time  nistories  agree  quite  well  with 
recorded  time  histories,  read  from  35  ram.  motion  picture  records 


* Numbers  refer  to  references  listed  at  end  of  report 
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A Linear  Theory  of  Ship  Motion  in  Irregular  Waves 


Introduction 

In  a report  by  Fuchs  and  Einarson^^j  the  motion  of  a freely  floating 
fora  in  periodic  progressive  waves  was  studied.  The  motions  of  heaving  and 
pitching  were  analysed  into  a series  of  powers  of  the  wave  steepness  and  the 
first  two  terms  were  given  explicitly  for  a rectangular  block.  The  prin- 
cipal effect  of  the  nonlinearitys  as  given  by  terns  of  the  order  of  the 
square  of  the  wave  steepness,  is  to  introduce  coupling  between  heaving, 
pitching  and  surging  and  to  introduce  small  terms  of  higher  frequencies, 
which  cause  a shift  of  maximum  values  and  an  additional  oscillation  in 
the  wave  trough. 

It  is  the  object  of  the  present  paper  to  extend  the  previous  theory 
so  that  it  applies  to  complex  wave  motions.  This  practically  restricts 
our  efforts  to  the  linear  theory,  and,  hence,  the  problem  is  a special 
case  of  the  linear  prediction  theory  discussed  by  Fucha(2).  Thus  the  oscil- 
lations of  a ship  in  a complex  system  of  progressive  waves  of  small  steep- 
nesses can  be  expressed  as  the  convolution  type  integral  of  the  recorded 
wave  notion  and  a kernel  function  which  is  the  Fourier  integral  of  the 
response  of  the  ship  to  a sinusoidal  forcing  function.  Kernels  for 
pitching  and  heaving  are  computed  explicitly  for  a freely  floating  rectan- 
gular block  using  experimental  values  of  added  masses  and  damping  coeffi- 
cients, as  determined  from  oscillations  in  still-  water.  These  kernel 
functions  are  applied  to  recorded  irregular  wave  motions  and  the  predictions 
are  compared  with  values  read  from  35  ram.  motion  picture  records  of  the 
motion.  The  experimental  work  has  been  reported  separately  by  Sifcul(3). 
Similar  investigations  have  been  made  with  a model  ship.  In  order  to 
avoid  the  laborious  numerical  integrations  involved  in  computing  the 
responae  functions,  these  were  determined  experimentally. 


Theory 

We  adopt  Kriloff's  original  auggestion  and  assume  that  the  force 
exerted  by  the  wavee  on  the  ehip'e  hull  can  be  represented  as  the  sum  of 
a virtual  mass  force,  a linear  damping  foree  and  the  force  arising  from 
the  incident  wave  acting  on  the  hull  in  its  equilibrium  position.  It 
is  assumed  that  the  virtual  mass  and  damping  coefficients  are  constants 
which  can  be  determined  from  oscillations  in  still,  water,  although  in 
reality  they  vary  with  frequency  as  has  been  shown  by  Havelock(4j,  with 
these  assumptions,  we  can  picture  the  ship  in  waves  as  a linear  dynamical 
aystea  for  which  the  forcing  function  is  the  sum  of  the  unperturbed  wave 
force  and  the  additional  force  due  to  oscillations  in  still  water,  with 
modifications  for  the  relative  motion.  Thus  we  explicitly  neglect  the 
effect  of  the  diffracted  wave.  Note  that  the  complete  solution  of  the 
linear  theory  problem  requiree  the  solution  of  the  diffraction  and  forced 
oscillations  problems. 
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2. 


Under  these  conditions  the  equations  of  hearing  and  pitching  are 
respectively 


~ (C  -n)dx+g/>A£  =jjp  cos  fly  dS 

s 


(1) 


+ p qB\?0  - jfjflxpcos  ny-y  pcosnx)dS  (2) 

The  second  and  third  terns  on  the  left  hand  sides  of  these  equations 
represent  the  vertical  nass  force  or  monent  of  inertia  and  the  damping 
forces,  respectively.  The  effect  of  relative  notion,  indicated  by  the 
difference  ( £ - n ) in  the  integrated  terms,  is  considered  only  for  heaving, 

since  the  additional  calculations  introduced  are  rather  tedious  and  do  not 
appreciably  effect  the  result.  Substituting  for  p the  unperturbed  wave 
pressure  we  show  in  Appendix  I that  the  heaving  and  pitching  motions  of  a 
rectangular  block  are  given  in  terms  of  the  wave  amplitude  17  by  complex 
responee  factors  r*  and  rp  in  the  forms 

£ * rh  (or ) 77  (<r  , t)  , 9 * rp  (<r)r)  {a,  t)  (3) 


The  functions  rh  (cr)  and  rp  (cr)  are  given  explicitly  on  pages  9 and  10, 

Appendix  II,  and  their  real  and  imaginary  parts  are  plotted  versus  cr 
in  Figures  1 and  2,  The  real  and  imaginary  parts  of  such  a complex 
response  factor  represent  "in-phase"  and  "out-phase"  components  respectively. 

For  a progressive  wave  motion  which  is  not  purely  sinusoidal,  the 
general  theory  presented  by  Puchs(2)  j.g  applicable  and  the  heaving  and 
pitching  motion  may  be  written  in  the  form: 


E(t ) 


J rj(T)  Kc  ( t -r  ) dr 


(4) 


where  the  kernel  is 


Kg  (t)  * J rg  {o-)e*lflrT  dr 


(5) 


rE  (<r)  being  the  reponse  of  the  block  to  a sinusoidal  wave  of  unit 
amplitude  and -frequency  cr  . Mote  that  r£  (cr)  has  real  and  imaginary 
parts  which  are  even  and  odd,  respectively,  which  guarantees  that  E is 
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real  for  a real  77  . Thue  if  we  write  rE  = /E  + i YE  , where  XE  , YE  are 
real  functions  of  a , then  KE  can  be  written  in  the  form 

KE(r)=^jfxt(a)  cos  order  YE  (r)  sin  <tt  dT  (6) 

= Ke(c3(  r ) + KEcs)(r) 

The  details  of  the  calculation  of  these  kernels  are  given  in  Appendix  II. 

The  results  are  plotted  in  Figures  3 and  4* 

The  principal  advantage  of  this  formulation  is  that  the  kernel  is 
independent  of  the  particular  wave  system  encountered  and  therefore  need 
be  calculated  only  once,  thereby  eliminating  repetitious  operations  other- 
wise required.  The  single  operation  of  Equation  (4)  combines  the  three 
separate  operations  of  the  common  Fourier  analysis  consisting  of  an  analysis 
of  v^Ct  ) into  its  Fourier  components,  multiplication  of  each  component  by 
its  response  factor  and  finally  synthesis  to  obtain  E(t). 

The  kernel  can  of  course  be  computed  if  the  response  is  determined 
empirically.  For  our  purpose  these  methods  are  convenient  for  ship  models, 
since  they  replace  the  laborious  calculation  of  force  components  due  to 
the  unperturbed  wave  pressure  acting  on  the  immersed  portion  of  the  hull, 
with  simple  laboratory  measurements.  This  calculation  can  be  performed 
either  by  a series  of  teste  of  the  model  in  sinusoidal  waves  or  by  Fourier 
analysis  of  records  secured  in  a wave  motion  having  a broad  freqnecy 
spectrum.  The  latter  method  has  been  carried  out  for  the  rectangular 
block  and  the  results  are  presented  in  Figure  5.  Run  8-A  was  used  to  com- 
pute the  heaving  response,  from  which  the  corresponding  heaving  kernel  is 
obtained  directly.  This  kernel  was  then  applied  to  Bun  8-A  as  a check, 
and  also  to  Run  10  which  did  not  enter  in  the  calculations  of  the  kernel. 

The  results  for  pitching  were  much  less  satisfactory,  since  the  im- 
portant range  of  frequneciea  in  the  pitching  response  occurs  for  higher 
frequencies  (see  Figures  1 and  2);  and  for  these  frequencies  the  trans- 
forms are  difficult  to  compute  accurately.  It  is  suggested  therefore  that 
this  technique  probably  could  be  used  to  advantage  only  for  heaving. 

It  should  be  noted,  finally,  that  the  calculations  will  be  more  accurate, 
the  more  irregular  the  wave  form  used  in  the  experiment;  since  only  for 
a fairly  irregular  wave  will  components  of  different  frequencies  be  present 
in  appreciable  amounts. 

The  second  method  consists  of  merely  placing  the  model  to  be  tested 
in  uniform  wave  trains  of  several  different  frequencies  and  recording  the 
heaving  and  pitching  motions.  Suppose  that  a surface,  profile  of  the 
forms 


7}  ~ d cos  er  t 


(7) 
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is  used  and  the  resulting  heaving  and  pitching  notions  are 
C - Co  cos  (O*  t - « ),  5 s <£0  cos  (<r  t - a) 


(8) 


then  Xh,  Yh,Xp,  Yp  may  be  determined  by  the  pairs  of  equations 

InJ  VXh2+Vh2  = £0/a  , V*h  = ten  * 

(9) 

|rPl=V/V"+Yp2  VXP  = r 

The  ambiguity  in  sign  may  be  eliminated  by  the  relations 


X h = I rh I cos  € , Yh  = I r I sin  € , Xp  = lrp  I cos  r,  (io) 

Yp  = I rp  I sin  t 

This  analysis  was  carried  out  on  an  actual  ship  model.  The  empirical 
reponse  factors  were  determined  and  are  plotted  in  Figures  10  and  11. 

Then  transforms  were  obtained  to  give  the  kernel  functions  for  heaving  and 
pitching.  The  kernels  were  then  applied  to  an  irregular  wave  train,  and 
the  results  are  shown  in  Figure  12. 

In  considering  the  effects  due  to  non-uniforra  wave  trains,'  a method 
which  immediately  suggests  itself  is  the  treatment  "wave  by  wave".  In 
this,  each  oscillation  of  the  surface  would  be  considered  as  part  of  an 
infinite  wave  train  of  height  and  period  equal  to  that  of  the  wave  in 
question.  The  results  of  this  report  offer  some  indications  of  the 
validity  of  this  approach.  To  this  end  such  an  analysis  was  carried 
out  on  two  of  the  runs,  namely  Huns  13  and  8-A.  The  results  are  indicated 
in  Figures  6 and  9 by  the  small  circles  which  represent  the  amplitude  and 
approximate  position  of  the  maximum  heaving  and  pitching.  It  is  seen 
that  while  the  results  are  fairly  good  for  Run  13,  which  is  relatively 
smooth,  this  analysis  does  not  give  even  a rough  approximation  in  Run  8-A, 
which  is  quite  irregular.  It  is  noted,  furthermore,  that  the  results  for 
pitching  are  poorer  than  the  results  for  heaving.  An  examination  of  the 
kernel  function  in  Figures  3 and  4 will  show  that  this  is  to  be  expected, 
for  the  pitching  kernel  dies  out  more  slowly  with  r . Generally  speaking 
the  rate  at  which  the  amplitude  of  the  kernel  of  a given  effect  decreases 
with  r is  a measure  of  the  dependence  of  that  effect,  at  a given  instant, 
on  the  past  and  future  behavior  of  the  surface  elevation. 

A second  problem  which  may  be  discussed  in  connection  with  this  report 
is  the  difference  in  the  central  frequency  of  the  Fourier  spectra  of  the 
surface  elevation  and  the  pitching.  If  one  has  a given  power  spectrum, 
that  is  the  sum  of  the  squares  of  the  amplitudes  of  the  cosine  and  sine 
Fourier  transforms,  of  a surface  elevation,  then  the  power  spectrum  for 
pitching  motion  may  be  obtained  from  it  by  multiplying  ty  the  absolute 
value  of  the  pitching  response  rp  (cr) . Since  this  latter  quantity  has 


its  maximum  at  a fixed  frequency,  and  the  surface  spectrum  may  be  cex:tered 
about  different  frequencies,  the  result  of  this  multiplication  would  be  to 
shift  the  central  frequency  toward  the  center  of  the  response  curve.;  that 
is,  a ship  would  tend  to  pitch  at  its  natural  frequency.  Two  example £ of 
this  process  are  presented  in  Figure  11,  u,-? mg  hypothetical  spectra  for  rj . 
The  absolute  value  of  the  pitching  response  is  plotted  in  Figure  13. 

These  examples  show  that  if  the  spectrum  of  tj  is  sufficiently  broad  the 

process  may  result  in  an  appreciable  shift  in  frequency.  The  effect  would 
also  be  present  in  heaving  but  would  probfebly  be  small  due  to  the  relatively 
flat  response  in  heaving.  This  appears  to  explain  the  observations  recently 
reported  by  Williams (o)  to  the  effect  that  under  confused  sea  conditions, 
the  predominant  rolling  period  was  the  natural  rolling  period.  It  does  not 
appear  to  suggest,  however,  the  absence  of  forced  rolling,  as  Haveloek(b) 
has  observed. 
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APPENDIX  I 


Calculation  of  Response  Factors  Tor  a Rectangular  Block 

In  a previous  report  by  Fuchs  and  Einarson(l),  the  response  factors 
were  computed  for  the  pitching  and  heaving  of  a freely  floating  rectangular 
block  in  periodic  waves.  However,  for  the  range  of  wave  lengths  considered 
in  this  report  a more  refined  analysis  of  the  relative  motion  terms  is 
required.  Owing  to  the  requirements  of  the  prediction  theory,  however, 
we  must  confine  our  attentions  to  the  linear  theory. 

The  forcing  functions  in  heave  and  pitch  are  then 


Fi 


A 

p cos  ny 


ds 


Lj 


xp  COS  ny  - yp  cos  (nx) ) ds 


(I.  1) 
(I. 


where  P is  the  pressure  field  of  the  incident  wave  and  the  integrations 
are  taken  over  the  immersed  portion  of  the  hull  when  it  is  in  its  still 
water  equilibrium  position.  For  the  range  of  frequencies  uith  whi.th  m 
shall  be  concerned,  there  is  no  appreciable  lack  of  accuracy  in  ascnu ming 

deep  water  conditions  for  all  components.  Then  we  have 


a 6 

P = p Q V-p-fiT 
where 

, go  ky  i(kx-crt) 

<f>  = i e e 


(I.  3) 


(I.  4) 


is  the  complex  velocity  potential  of  the  incident  waves  associated  with 
the  surface  elevation 


(I.  5) 


This  is  a periodic  wave  of  amplitude  a with  its  frequency  <7  and  wave 
number  k related  by  the  equation 

o-2  = gk  (I.  6) 

This  wave  moves  in  the  positive  x direction  with  the  constant  velocity 

c = <r/k  = Vg/k  •, 
for  convenience  we  write 

F,  - J p cos  ny  ds  - J p ( s)  d A0 

'S  *• 


(I.  7) 
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7. 


where  A0  is  the  projection  of  the  immersed  hull  surface  S0  on  the  xz 
plane.  Substituting,  we  find 


F|  zj  [g/>  (d+£+  B%)  + p ga  e kd  e'^k*  j d Ao 

Ae 


-kd  sin  k L -tat 
e 


(I.  8) 


: P Q V0  + />  g A0£  ~p  gBae'  k 

where  V0  is  the  immersed  volume  of  the  block  in  equilibrium,  A0  is  vh* 

equilibrium  value  of  the  water  plane  area,  and  B , L , d are  the  beam, 
half-length  and  draft,  respectively. 


Lj  = / x p d A o - / ypdP0 

« A-  J P- 


/A.  •'P. 

where  P0  is  the  projection  of  S0  onto  the  yz  -plane.  Substituting 
we  find 


(I.  9) 


Lf  - f \qp  ( d+£  + 9% ) + p q a e*kd  e'^x‘trt)j  ( xd  A0  - y d P0) 

- | pgBL5  ^-2iB^>gaek4d  [sin  kL-kLcoskL  e"at 


+ 2iB/>ga  J£rLk  (I  - h d e "a  - e*"0 ) e 
The  relative  motion  terms  are 
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~L  d' 


V A 
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d2  £ . z 2 sin  kL  - fcr t 

= M — r + M acr  — e 


(I.  10) 


(I.  11) 


d C , • 2 s’m  kL  -iat 

= Ni  + •O’O j, e 


(I.  12) 
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On  introducing  the  notation 
2b=  N|/M  , u)h2  -- 


2J.  P l 
M 


/•>'  - M/ML 


....  o m /-r  2 2 g B L3  2B  P g 

y 4, /ML,  2 c = Nz / 1 , wp2  = - - * x , /x  - 

p(a ) = - d-kk°  {rin  kL  - kL  cos  kL^  + 


{ i - kd  e K"  - e M] 


the  equations  can  be  written  in  the  form 

e 4.  i. 

: £ = | a h(cr) 

6 - i /i  p (cr)  ie"ltr| 


(I.  13) 


(I.  14) 
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d2£> 

dB 
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dt 

+ OJ  p 

n 2 , sm  kL  "I  cr  t 

e - (3cr  - i - -k--  /cr  | ae 


-kd  Sir  k L 
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(I.  15) 


neglecting  the  transient  oscillation  which  dies  out  exponentially,  the 
forced  oscillation  can  be  written  in  the  form 


£ (t  ) = r h (o')  Tj  ( t ) 

P (t  ) = rp  (cr)  7 7 (cr) 


Cl.  16) 

(I.  17) 


where  r h and  rp  are  complex  response  factors  given  in  Appendix  II. 

These  functions  are  plotted  in  Figures  1 and  2. 

The  transition  to  irregular  wave  forms  is  accomplished  by  respresenting 
the  surface  by  its  Fourier  transform 

, +at  __ 

(I.  18) 


rj{  \ ) - I /*/2 7T  / a (cr)  e’ 1 ^ 1 d cr 

t'—  m 


The  preceding  results  may  then  be  applied  by  making  use  of  the  principle 
of  superposition  to  give,  for  heaving  or  pitching. 


CD 

rE  (cr)  a (cr)  e"llTtd  cr 
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(I.  19) 


This  last  integral  may  then  be  transformed  into  a ccr.f  elution  integral 
involving  the  surface  elevation  and  a kernel 

,+0D 
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where 


/’’■CD 

77(r)KE(t-r)dT 


Kp (r)  = I 


/-/Ztt  [ 

J-QC. 


rF  (cr)  e"  crr  d cr 


(I.  20) 


(I.  21) 


It  has  already  been  indicated  that  the  response  may  be  determined 
empirically  by  measurements  in  an  Irregular  wave  train.  To  see  this  we 
take  the  Fourier  transform  of  equation  (I.  19).  Then  noting  tnat  th« 
response  rE  is  the  inverse  transform  of  the  kernel  Kg  , we  obtain 


rE  (it) 


bg  (cr) 
o (<r) 


(I.  22^ 


where  bg(cr)  and  a (cr)  are  the  transforms  of  the  heaving  or  pitching 
and  surface  elevation  respectively. 


APPENDIX  II 


Numerical  Calculations 


For  convenience,  the  equations  used  to  determine  the  heaving  and 
pitching  kernels  together  with  the  values  of  the  numerical  constants  in- 
volved are  list*'  below,  X and  Y representing  the  real  and  imaginary 
parts  of  r , respectively. 


sQD  * - s© 

Kh  (r)  = v/2 Ttt J Xh(cr ) cos  cr  r d<r  + J Yh  (cr)  sin  err  do-  =Khc(r  )+Kh*(r) 
Xh  (<r) 


(ah((T).  mkL  £Cr*){Cj*-(?)  4-  2b  a* 


(wh2-o-2)2+  4b2  <tz 


(II.  1) 


sin  kl. 


Yh  (O’)  = 


2b  ><r  n(o-) -Zbfi  a3  - y k 


(<oh2-  a2) 


(wh"  - o-2)2+4  b‘  <t£ 


r h = ‘'/^Xh^~+  Y h 2 
h(o-)  * ekd 

L = 0.50  ft.,  d = 0.109  ft.,  B = 0.425  ft.,  M = 0.199  lbs.  ft”1  sec.2, 

M = 0,097  lbs.  ft. “I  sec.2,  N|  = 0.66  lbs.  ft.-1  sec.,  Wh2  = 133  sec.”2, 

2Brg 

2b=N|/M=  3,30  sec."1,  a = = 267  ft."1  sec.”2,  f3  s M/ML  = 1.00  ft.-1, 

j-Ni/ML=6.66  ft.”1  sec.”1. 


sec 


10 


CL 

cr 

= y 2/rr  f Xp(o-)  cos  order  ■ 
Jo 

h-s/2 lir  Y 

p(cr)sin<7Tdcr  = Kpc(r)+ 

Kp*(r) 

Xp  (cr) 

~2c  cr  p (cr) 

(cjp  2 - cr 2 )z  + 4cz  a-2 

(II. 

2) 

Yp(cr) 

fi  pier) 

(GDp2-^*)*  + 4c2  cr2 

r p : ■v/  XpZ  + Yp  2 

p(o-)  = 

k2  |sin  kL  - kL  cos  kl_j 

sin  k L 
k2 

| ! - k d e‘kd-e'kd  j 

I = 0.014 

lbs.  ft.  sec. Nz= 0.071  lbs.  ft. 

sec.,  ^pz  = 133  sec.' 

-2 

J 

P I 

3780  ft. ”3  Sec.“2,  2c  = 

Ng/I =5.26 

sec. 

Two  methods  of  computation  were  used  to  carry  out  the  integration  in 
equations  (II.  1)  and  (II.  2).  In  keeping  with  the  approximate  nature  of 
the  pitching  response  due  to  neglect  of  the  relative  motion  terms,  an 
approximate  method  was  used  for  equation  (II.  2).  If  the  function  Xp (<r) 

is  extended  so  as  to  form  an  even  function  of  cr  , then  over  the  range 

-a  <a  <+  cr  it  may  be  approximated  bj  a finite  Fourier  sum  involving 
onlj  cosine  termsj  that  is, 

Xp  (cr)  = a0  + a,  cos  = - + . . . + an  cos  —^5  (II.  3) 

where 


Ok  : 


2_ 

a 


n TTCT 
W 


d cr 


(n.  4) 


Hence,  if  cr  is  taken  sufficiently  large  that  Xp(cr)  is  approximately 
zero  for  cr>o,  the  Ok's  give  an  approximation  to  Kp (cr)  at  the  points 

k7T  , 

0-  , k 0)..<10  j that  3.8, 


(II.  5) 


In  an  exactly  analogous  manner,  if  Yp(<r ) is  extended  so  as  to  form  an 
even  function  it  may  be  approximated  by  a finite  sine  series,  with  co- 
efficients b| , . . . , bn  , and  these  give 


bk 


k = I... . ,n 


(II.  6) 


~7 


11. 


A Fourier  analysis  was  carried  out  using  the  standard  24  ordinate  scheme, 
giving  the  value  of  Kpc  at  thirteen  values  of  <r  , and  Kp*  at  eleven 
values. 

The  integrals  in  equation  (II.  1)  were  evaluated  with  somewhat  more 
accuracy.  For  large  values  of  t asymptotic  representations  are  available 
for  KhC  (r)  and  K^Cr)  in  the  form 


,,  c / _ \ . Xh*  (0)  . X'"  (o) 

Kh  (r ) = - — —1 + — r 


(0) 


+ - 


, Y„(o)  Yh"(o)  Yh'V  (o) 

Kh  (r)  - — - - -TJ  - ~9- — - . . . 


(u.  7) 


where  primes  indicate  differentiation  with  respect  to  O’  . The  indicated 
differentiations  were  carried  out  and  it  was  found  that  K*,c  and  Khs 

may  be  taken  to  be  zero  for  t >0.8.  In  the  region  t<0.8,  Khc  and 

Khs  were  computed  using  Filon's  formula^)  for  trigonometric  integrals. 

The  errors  in  the  use  of  this  formula  can  be  estimated  from  the  third 
differences  of  X^o-)  and  Y^o-),  and  were  kept  less  than  two  per  cent 

of  the  maximum  values. 


The  convolution  integrals  involving  7)  ( t ) and  tne  kernels  for  neaving 
and  pitching  were  carried  out  using  the  trapezoid  rule  with  an  interval  of 
0.1  seconds.  The  integrations  involved  in  determining  the  transforms  of  f) 
and  £ , which  led  to  the  empirical  response  factor  for  heaving  were  carried 
out  using  Filon's  formula  and  the  kernel  obtained  by  applying  the  Fourier 
series  technique  already  outlined. 
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FIGURE  I 


COMPONENTS  OF  THE  HEAVING  RESPONSE  OF  A RECTANGULAR  BLOCK 


HEAVING  MOTION  USING  AN  EMPIRICALLY  DETERMINED  HEAVING  KERNEL 
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PITCHING  AND  HEAVING  OF  A RECTANGULAR  BLOCK  IN  A WAVE  GROUP 


analysis  of  pitching  and  heaving  of  a rectangular  block  in  irregular  waves 


ANALYSIS  OF  PITCHING  AND  HEAVING  OF  A RECTANGULAR  BLOCK  IN  IRREGULAR  WAVES 
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EMPIRICALLY  DETERMINED  RESPONSE  FACTORS  FOR  HEAVING  OF  SHIP  MODEL 


